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, Abstract 

The overwhelming majority of survivable (fault-tolerant) network design models assume a 
uniform scenario set. Such a scenario set assumes that every subset of the network resources 
^vq ■ (edges or vertices) of a given cardinality k comprises a scenario. While this approach yields 

problems with clean combinatorial structure and good algorithms, it often fails to capture the 
true nature of the scenario set coming from applications. 

One natural refinement of the uniform model is obtained by partitioning the set of resources 
into faulty and secure resources. The scenario set contains every subset of at most k faulty 
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^ resources. This work studies the Fault- Tolerant Path (FTP) problem, the counterpart of the 

Shortest Path problem in this failure model. We present complexity results alongside exact and 
approximation algorithms for FTP. We emphasize the vast increase in the complexity of the 
problem with respect to its uniform analogue, the Edge-Disjoint Paths problem. 

^ ■ 1 Introduction 

\D 

The Minimum-Cost Edge-Disjoint Path (EDP) problem is a classical network design problem, 
defined as follows. Given an edge- weighted graph G = (V, E), two terminals s, t € V and an integer 
parameter k € Z+, find k edge-disjoint paths connecting s and t with minimum total cost. EDP is 
motivated by the following survivable network design problem: what is the connection cost of two 
nodes in a network, given that any k — 1 edges can be a-posteriori removed from the graph. The 
implicit assumption in EDP is that every edge in the graph in equally vulnerable. Unfortunately, 
this assumption is unrealistic in many applications, hence resulting in overly-conservative solutions. 
This paper studies a natural refinement of the EDP problem called the Fault- Tolerant Path Problem 
(FTP), formally defined as follows. 



Instance: An edge-weighted graph G = (V,E), two nodes s,t € V, a subset of the edges 
M C E and an integer k G Z + . 

Problem: Find minimum cost S C E such that S\F contains an s-t path for every 
F C M with \F\ < k. 



Note that FTP becomes EDP when M = E. We let fi(M, k) denote the set of possible faults, 
namely fi(M, k) = {A C M : \A\ < k}. 

A well-known polynomial algorithm for EDP works as follows. Assign unit capacities to all 
edges in G and find a minimum-cost A;- flow from s to t. The integrality property of the Minimum- 
Cost Flow (MCF) problem guarantees that an extreme-point optimal solution is integral, hence it 
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corresponds to a subset of edges. It is then straightforward to verify that this set is an optimal 
solution of the EDP problem (for a thorough treatment of this method we refer to the book of 
Schrijver [19]). 

The latter algorithm raises two immediate questions concerning FTP. The first question is 
whether FTP admits a polynomial time algorithm. In this paper we give a negative answer to this 
question conditioned on P^NP, showing that FTP is NP-hard. In fact, the existence of constant- 
factor approximation algorithms is unlikely already for the restricted case of directed acyclic graphs. 
Consequently, it is natural to ask whether polynomial algorithms can be obtained for restricted 
variants of FTP. We devote the first part of this paper to this question. 

The second question concerns the natural fractional relaxation of FTP. As we previously ob- 
served, one natural relaxation of EDP is the MCF problem. This relaxation admits an integrality 
gap of one, namely the optimal integral solution value is always equal to the corresponding optimal 
fractional value. In the context of FTP it is natural to study the following fractional relaxation, 
which we denote by Fractional FTP (FRAC-FTP). 



Instance: An edge weighted graph G = (V, E), two nodes s, t G V, a subset of the edges 
M C E and an integer k G Z + . 

Problem: Find a minimum cost capacity vector x : E — >■ [0, 1] such that for every 
F G Q(M, k), the maximum s-t flow in Gf = (V, E \ F), capacitated by x, is at least one. 



It is obvious that by further requiring x G {0, 1} E we obtain the FTP problem, thus FRAC-FTP 
is indeed a fractional relaxation of FTP. We later show that the integrality gap with respect to 
this relaxation is bounded by k + 1. Furthermore, we show that this bound is tight, namely that 
there exists an infinite family of instances with integrality gap arbitrarily close to k + 1. This result 
also leads to a simple (k + l)-approximation algorithm for FTP, which we later combine with an 
algorithm for the case k = 1 to obtain a ^-approximation algorithm. 

The remainder of the paper is organized as follows. In the following section we review some 
related work. In Section [3] we study the complexity of FTP. In Section [4] we provide some exact 
polynomial algorithms for special cases of FTP. In Section [5] we relate FTP and FRAC-FTP by 
proving a tight bound on the interagrality gap and show how this result leads to a /c-approximation 
algorithm for FTP. In Section [6] we conclude the paper and mention some interesting directions for 
future work. 

2 Related Work 

The shortest path problem was studied in numerous robust settings ( [201 EEJ IS E21 [21 [HI [6] ) . Yu 
and Yang [20] studied the shortest path problem with a discrete scenario set in the minimax and 
minimax-regret models. They show that both problems are NP-Hard even for 2 scenarios and 
layered graphs. Aissi et al. [2] extended these results by showing that both version of the problem 
admit an FPTAS if the number of scenarios is bounded. They also show that both problems 
are not approximable within 2 — e for any e > with an unbounded number of scenarios unless 
P=NP. Adjiashvili et al. [1] studied a two-stage feasibility problem in which for a given pair of 
terminals s and f in a graph G = (V,E), and two parameters /, r G Z+, the goal is to choose 
a minimum-cardinality subset of the edges S C E with the property that for every failure set 
F C E with \F\ < f there exists a recovery set R C E \ F with \R\ < r, such that s and t are 
connected in (S \ F) U R. The authors give an exact algorithm for f = r = 1 and a 2-approximation 
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algorithm for / = 1 and arbitrary r. Also some complexity results are provided. Puhl [18] studied 
the recoverable-robust shortest path problem, a robust path problem in a two-stage optimization 
setting. Hardness results were provided for discrete scenario sets, interval scenarios and T-scenarios 
for two variants of cost function. For one variant an approximation algorithm is given. Busing [6 
developed approximation algorithms and proved hardness results for the problem of finding the 
smallest number of edges in the graph which contains a shortest path according to every scenario. 

The robustness model proposed in this paper is natural for various classical combinatorial 
optimization problems. Of particular interest is the counterpart of the Minimum Spanning Tree 
problem. This problem is closely related to the Minimum A;-Edge Connected Spanning Subgraph 
(ECSS) problem, a well-understood robust connection problem. Gabow, Goemans, Tardos and 
Williamson [llj developed a polynomial time (1 + ^-approximation algorithm for ECSS, for some 
fixed constant c. The authors also show that for some constant d < c, the existence of a polynomial 
time (1 + ^-approximation algorithm implies P=NP. An intriguing property of ECSS is that the 
problem becomes easier to approximate when k grows. Concretely, while for every fixed k, ECSS is 
NP-hard to approximate within some factor aj. > 1, the latter result asserts that there is function 
j3{k) tending to one as k tends to infinity such that ECSS is approximable within a factor j3{k). 
This phenomenon was already discovered by Cheriyan and Thurimella [8], who gave algorithms 
with a weaker approximation guarantee. The more general Generalized Steiner Network problem 
admits a polynomial 2-approximation algorithm due to Jain [14J. 

We review next some of the related research in robust combinatorial optimization. Liebchen 
et al. [IT] proposed the notion of recoverable robustness as a tractable framework for robustness in 
railway optimization. The authors present a general recipe for obtaining two-stage robust problems 
from any combinatorial optimization problem. In [3] Ben-Tal et al. proposed a method to delay 
some decisions in an uncertain convex optimization problem to the stage at which the real data 
is revealed. Their approach does not force some decision variables X to be fixed as first-stage 
decision variables. Instead, these variables are allowed to be associated in the first stage decision 
with affine functions A"(C) of the uncertain input vector The resulting formulation turns out 
to be theoretically and practically tractable in many cases. The framework of Bertsimas and 
Caramanis [4] allows to deal with discrete changes in the second stage decision in the special case 
of linear optimization. 

Kouvalis and Yu [16j provide an extensive study of complexity and algorithms for robust coun- 
terparts of many classical discrete optimization. Instead of resource removal, these models assume 
uncertainty in final cost of the resources. Bertsimas and Sim [5] proposed a framework for incor- 
porating robustness in discrete optimization. They show that robust counterparts of mixed integer 
linear programs with uncertainty in the matrix coefficients can be solved using moderately larger 
mixed integer linear programs with no uncertainty. 

A two-stage robust model c&lleddemand-robust optimization was introduced by Dhamdhere et 
al. [9]. In this model a scenario corresponds to a subset of the initially given constraints that need to 
be satisfied. For shortest path this means that the source node and the target node are not fixed, 
but are rather revealed in the second stage. The authors provide approximation algorithms for 
some problems such as Steiner Tree, Multi-Cut, facility Location etc. in the demand-robust setting. 
Golovin et al. |12] later improved some of those results, including a constant factor approximation 
algorithm for the robust shortest path problem. Feige et al. [10] and Khandekar et al. [15] studied 
an extension of the demand-robust model that allows exponentially many scenarios. 
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3 Complexity of FTP 



Our first observation is that FTP generalizes the Directed m-Steiner Tree Problem (m-DST). The 
input to m-DST is a weighted directed graph G = (V,E), a source node s G V, a collections of 
terminals T C V and an integer m < \T\. The goal is to find a minimum-cost arboresence A C E 
rooted at s, that contains a directed path from s to some subset of m terminal. 

The m-DST is seen to be a special case of FTP as follows. Given an instance / = (G, s, T, m) 
of m-DST define the following instance of FTP. The graph G is augmented by \T\ new zero-cost 
directed edges E' connecting every terminal u 6 T to a new node t. Finally, we set M = E' and 
k = m — 1. The goal is to find a fault-tolerant path from s to t in the new graph. It is now 
straightforward to see that a solution S to the FTP instance is feasible if and only if S n E contains 
a feasible solution to the /c-DST problem (we can assume that all edges in E' are in any solution 
to the FTP instance). 

The latter observation implies an immediate conditional lower bound on the approximability of 
FTP. Halperin and Krauthgamer [13] showed that m-DST cannot be approximated within a factor 
log 2_e m for every e > 0, unless NP C ZTlME(n polyl ° 9 ^). As a result we obtain the following. 

Proposition 1. For every e > FTP admits no polynomial approximation algorithms with ratio 
log 2 ~ e k, unless NP C ZTIME(n polyl ° 9 ^ ). 

Remark 2. The reduction above can be easily adapted to obtain a ^-approximation algorithm for 
FTP for the special case that M C {e 6 E : t E e} using the algorithm of Charikar et. al. [TJ. 

In fact, any approximation algorithm with factor p(k) for FTP is an approximation algorithm 
with factor p(m) for m-DST. The best known algorithm for the latter problem is due to Charikar 
et. al. [7J. Their result is an approximation scheme attaining the approximation factor of m e for 
every e > 0. 

In light of the previous discussion, it seems that obtaining an approximation algorithm for FTP 
with a poly-logarithmic approximation guarantee is a challenging task. In fact, although we have 
no concrete result to support this possibility, such an algorithm might not exist. We end this 
discussion with another special case of FTP, which we denote by Simultaneous Directed m-Steiner 
Tree (m-SDST). An input to m-SDST specifies two edge-weighted graphs G\ = (V,Ei,wi) and 
G2 = (V, E2,W2) on the same set of vertices V, a source s, a set of terminals T C V and an integer 
m < \T\. The goal is to find a subset U C T of m terminals and two arboresences S\ C E\ and 
£2 C E2 connecting s to U in the respective graphs, so as to minimize w\(S\) + 102(82)- m-SDST 
is seen to be a special case of FTP via the following reduction. Given an instance of m-SDST, 
construct a graph G = (V',E) as follows. Take a disjoint union of G\ and G2, where the direction 
of every edge in G2 is reversed. Connect every copy of a terminal u £ T in G\ to its corresponding 
copy in G2 with an additional zero-cost edge e u . Finally, set M = {e u : u € T} and k = m — 1. A 
fault-tolerant path connecting the copy of s in G\ to the copy of s in G2 corresponds to a feasible 
solution to the m-SDST instance with the same cost, and vise- verse. 

4 Polynomial Special Cases 

This section is concerned with tractable restrictions of FTP. Concretely we give polynomial al- 
gorithms for arbitrary graphs and k = 1, directed acyclic graphs (DAGs) and fixed k and for 
Series-parallel graphs. We denote the problem FTP restricted to instances with some fixed k by 
fc-FTP. 
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4.1 1-FTP 



For simplicity of the exposition we assume that G is undirected. For the adaptation to directed 
graphs see Remark [71 We start with the following useful structural lemma. 

Lemma 3. Let X* be an optimal solution to FTP on the instance (G,M, k). Let H' = (V',X') be 
the graph obtained by contracting all edges of X* \M in H* = (V,X*). Let s' and t' be the vertices 
in V' , to which s and t were contracted, respectively. Then either s' = t' , or there are k + 1 edge 
disjoint s'-t' paths in H' . 

Proof. If X* C E\M then s' = t' . Assume this is not the case, namely that X* n M / 0. Observe 
that H' must contain at least k + 1 edge disjoint s'-t' paths: indeed if there was a cut of cardinality 
k or less in H, this would also be a cut in H*. Since X' C M, it corresponds to a failure scenario, 
namely X' € Q(M,k). □ 

Definition 4. An s-t bipath in the graph G = (V, E) is a union of two s-t paths Pi, P2 C E, such 
that for every two nodes u, v incident to both P\ and P2, the order in which they appear on Pi and 
P2, respectively, when traversed from s to t is the same. 

In the context of 1-FTP, we call a bipath Q = P 1 U P 2 robust, if it holds that Pi n P 2 D M = 0. 
Note that every robust s-t bipath Q in G is a feasible solution to the 1-FTP instance. Indeed, 
consider any failure edge e £ M. Since e Pi n P2 it holds that either Pi C Q — e, or P2 C Q — e. 
It follows that Q — e contains some s-t path. The next lemma shows that every feasible solution of 
the 1-FTP instance contains a robust s-t bipath. 

Lemma 5. Every feasible solution S* to the 1-FTP contains a robust s-t bipath. 

Proof. We assume without loss of generality that S* is a minimal feasible solution with respect to 
inclusion. Let Y C S* be the set of bridges in (V, S*). From feasibility of S* , we have Y n M = 0. 
Consider any s-t path P in S* . Let u\, ■ ■ ■ , u r be be the set of vertices incident to Y = Pfl Y. Let 
U{ and Ui+i be such that Uitii+i Y. (if such an edge does not exist, we have Y = P, which means 
that P is a robust s-t bipath). Note that S* must contain two edge-disjoint Ui-u; L+ i paths L\,L2. 
Taking as the set Y together with all such pairs of paths L±,L2 results in a robust bipath. □ 

We can conclude from the previous discussion and Lemma [5] that all minimal feasible solutions 
to the 1-FTP instance are robust bipaths. This observations leads to the simple algorithm, which 
is formally given as Algorithm [TJ We prove its correctness in the following theorem. 

Theorem 6. There is a polynomial algorithm for 1-FTP. 

Proof. To solve 1-FTP we need to find the minimum cost robust s-t bipath. To this end let us define 
two length functions £±,£2 ■ V x V — > M+. For two vertices u, v € V let £i(u, v) denote the shortest 
path distance from u to v in the graph (V, E\M), and let £2{u, v) denote the cost of the shortest pair 
of edge-disjoint u-v paths in G. Clearly, both length functions can be computed in polynomial time 
(e.g. using flow techniques). Finally, set £(u,v) = mm{£i(u,v),£2(u,v)}. Construct the complete 
graph on the vertex set V and associate the length function £ with it. Observe that by definition 
of £, any s-t path in this graph corresponds to a robust s-t bipath with the same cost, and vice 
versa. It remains to find the shortest s-t bipath by performing a single shortest s-t path in the new 
graph. For every edge uv in this shortest path, the optimal bipath contains the shortest u-v path 
in (V, E \ M) if £(u, v) = £i(u, v), and the shortest pair of u-v paths in G, otherwise. □ 
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Algorithm 1 : FindBipath(G, s, t, M) 



Input: G = (V,E) a directed graph, s,t G V and M C E. 
Output: A minimum-cost s-t bipath Q. 
1: Q^0. 

2: G S AFE^(V,E\M). 
3: for all it, v G V x V do 

4: Compute £i(u, u), the shortest u-v path distance in Gsafe- 

5: Compute l2(u,v), the cost of the minimum u-u 2- flow in G with unit capacities. 

6: £(u,v) <— mm{£i(u,v),£2(u,v)}. 

7: Compute the shortest s-t path P in (1/, V x V) and length function t. 

8: for all uv G P do 

9: if £(u, u) = ^i(n, v) then 
10: Add to Q a shortest u-u path in Gsafe- 
11: else 

12: Add to Q a shortest edge-disjoint pair of u-v paths in G. 
13: return Q. 



Remark 7. Let us remark that Algorithm Q] easily extends to work on directed graphs as follows. 
Lemma [3] is replaced with the more general statement that the minimum s-t flow in the graph 
(V, X*) capacitated by the vector c e = 1 if e G M and c e = oo, otherwise is at least k + 1. The rest 
of the analysis remains the same. 



4.2 /c-FTP and Directed Acyclic Graphs 

Let us first consider the case of a layered graph. The generalization to a directed acyclic graph 
is done via a standard transformation, which we describe later. Recall that a layered graph G = 
(V, E) is a graph with a partitioned vertex set V = V\ U • - • U V r and a set of edges satisfying 
E C Uie[r-i] x Vi+i- We assume without loss of generality that V\ = {s} and V r = {t}. For 
every i G [r — 1] we let Ei = E n x T^+x- 

Analogously to the algorithm in the previous section, we reduce fc-FTP to a shortest path 
problem in a larger graph. The following definition sets the stage for the algorithm. 

Definition 8. An i- configuration is a vector d € {0, 1, • • • , k + l} Vi satisfying X^ey d v = k + 1. 
We let supp{d) = {v G V\ : d v > 0}. For an i-configuration d 1 and an (i + l)-configuration d 2 we let 

^(d 1 ,^ 2 ) = supp{d l ) Usupp(d 2 ) and E(d 1 ,d 2 ) = E[V(d l ,d 2 )]. 

We say that an ^-configuration d 1 precedes an {i + l)-configuration <i 2 if the following flow 
problem is feasible. The graph is defined as H(d ,d 2 ) = (V(d , d 2 ), E(d , d 2 )). The demand vector 
v and the capacity vector c are given by 

\-d\ if u G suppid 1 ) , f 1 if e G M 

= < o and c e = < 

ld„ if ii G supp(d l ) yoo if eeE\M, 

respectively. If d 1 precedes <i 2 we say that the link (d l ,d 2 ) exists. Finally, the cost £(d 1 ,d 2 ) of this 
link is set to be minimum value w(E') over all E' C E(d 1 ,d 2 ), for which the previous flow problem 
is feasible, when restricted to the set of edges E' . 
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The algorithm constructs a layered graph % = (V, 6) with r layers Vi, • • • ,V r . For every i £ [r] 
the set of vertices Vi contains all i-configurations. Observe that Vi and V r contain one vertex each, 
denoted by c s and c*, respectively. The edges correspond to links between configurations. Every 
edge is directed from the configuration with the lower index to the one with the higher index. The 
cost is set according to Definition [8j The following lemma provides the required observation, which 
immediately leads to a polynomial algorithm. 

Lemma 9. Every c s -c l path P in H corresponds to a fault-tolerant path S with w(S) < £{P), and 
vise-versa. 

Proof. Consider first a fault-tolerant path S C E. We construct a corresponding c s -c* path in H 
as follows. Consider any k + 1 s-t flow f s , induced by S. Let p l , ■ ■ ■ ,p l be a path decomposition 
of f s and let 1 < pi, ■ ■ ■ , p\ < k + 1 (with J2ie[i] Pi = k + 1) be the corresponding flow values. 

Since G is layered, the path p> contains exactly one vertex v\ from Vi and one edge from E% 
for every j € [/] and i € [r]. For every i S [r] define the i-configuration d l with 

4 = Yl pu 

if some path p 3 contains v, and d l v = 0, otherwise. The fact that d l is an z-configuration follows 
immediately from the fact that f s is a (k + l)-flow. In addition, for the same reason d % precedes 
di + i for every i S [r — 1]. From the latter observations and the fact that d 1 = c s and d r = c* it 
follows that P = d 1 , d 2 , ■ ■ ■ , d T is a c s — c* path in if with cost £(P) < w(S). 

Consider next an c s - c* path P = d 1 ,--- ,d r with cost £(P) = TJlZl ^{d\ d i+1 ) . The cost 
£(d l ,d l+1 ) is realized by some set of edges Ri C E(d l ,d' l+1 ) for every i € [r — 1]. From Definition [HI 
the maximal s-t flow in the graph G' = (V,R) is at least k + 1, where -R = Uj e r r _u-Ri- Next, 
Lemma [3J guarantees that there exists some feasible solution S C R, the cost of which is at most 
£(P). In the latter claim we used the disjointness of the sets Ri, which is due the layered structure 
of the graph G. This concludes the proof of the lemma. □ 

Finally, we observe that the number of configurations is bounded by 0(n k+1 ), which implies 
that /c-FTP can be solved in polynomial time on layered graphs. 

To obtain the same result for directed acyclic graphs we perform the following transformation of 
the graph. Let v\, ■ ■ ■ , v n be a topological sorting of the vertices in G. Replace every edge e = ViVj 
(i < j) with a path p e = Vi,uf +1 , ■ ■ ■ , Uj_ 1 ,Vj of length j — i + 1 by subdividing it sufficiently many 
times. Set the cost of the first edge on the path to w'(viuf +1 ) = w{viVj) and set the costs of all 
other edges on the path to zero. In addition, create a new set of faulty edges M' , which contains 
all edges in a path p e if e € M. 

It is straightforward to see that the new instance of FTP is equivalent to the original one, while 
the obtained graph after the transformation is layered. We summarize with the following theorem. 

Theorem 10. There is a polynomial algorithm for k-FTP restricted to instances with a directed 
acyclic graph. 

4.3 Series-Parallel Graphs 

Recall that a graph is called series-parallel (SRP) with terminal s and t if it can be composed from 
a collection of disjoint edges using the series and parallel compositions. The series composition 
of two SRP graphs with terminals s, t and s',t' respectively, takes the disjoint union of the two 
graphs, and identifies t with s'. The parallel composition takes the disjoint union of the two graphs 
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and identifies s with s' and t with t'. Given a SRP graph it is easy to obtain the aforementioned 
decomposition. 

The algorithm we present has linear running time whenever the robustness parameter k is fixed. 
The algorithm is given as Algorithm [2l In fact, the algorithms computes the optimal solutions SV 
for all parameters < k! < k. The symbol _L is returned if the problem is infeasible. 

Algorithm 2 : FTP-SeriesParallel(G, s, t, k) 

Input: G = (V,E) a series-parallel graph, s,t € V and M C E, k € Z + . 
Output: Optimal solution to FTP for parameters 0, 1, • • • ,k. 

1: if E = {e} A e G M then 

2: Return ({e}, , _L) 

3: if £ = {e} A e £ M then 

4: Return ({e}, • • • , {e}) 

=>■ G is a composition of H\,H2- 

5: (S(J, • • • , Si) <- FTP-SeriesParallel(#i, M n fc) 

6: (Sg, • • • , Sf) <- FTP-SeriesParallel(# 2 , M n £?[fl" 2 ], *) 

7: if G is a series composition of H\, H% then 

8: for i = 0, ■ ■ ■ , A; do 

9: if Sj =± V Sf =± then 

10: S{ 4-± 

11: else 

12: Si <- Sj U S? 

13: if G is a parallel composition of H\,H2 then 

14: mi max{i : S] ^-L} 

15: 771,2 ^— max{« : S'j 2 ^-L} 

16: for i = 0, ■ ■ ■ , A; do 

17: if i > mi + 7712 + 1 then 

18: S{ <-± 

19: else 

20: r <- argmin-xKjKiiwiSj) + wiSf^^)} // Sl x = S 2 _ x := 

21: S{ ^ 5*^ U S r ^_ r ^ 

22: Return (5q, • • • , 5^) 



Theorem 11. Algorithm® returns an optimal solution to the FTP problem on SRP graphs. The 
running time of Algorithm® is 0(nk). 

Proof. The proof of correctness is by induction on the depth of the recursion in Algorithm [2j 
Clearly the result returned by Algorithm [2] in lines 1-6 is optimal. Assume next that the algorithm 
computed correctly all optimal solutions for the subgraphs Hi , H2 , namely that for every i € [2] 
and j £ [k], the set Sj computed in lines 7-8 is an optimal solution to the problem on instance 
if = (H i ,MnE[H i ],j). 

Assume first that G is a series composition of Hi and H2, and let < i < k. If either Sj =_L or 
Sf =_L the problem with parameter i is clearly also infeasible, hence the algorithm works correctly 
in this case. Furthermore, since G contains a cut vertex (the terminal node, which is in common to 
Hi and H2), a solution S to the problem is feasible for G if and only if it is a union of two feasible 
solutions for Hi and H2- From the inductive hypothesis it follows that Si is computed correctly in 
line 14. 
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Assume next that G is a series composition of H\ and H 2 . Consider any feasible solution S' to 
the problem on G with parameter i. Let S[ and S 2 be the restrictions of S' to edges of H\ and 
H 2 respectively, and let n\ and 112 be the maximal integers such that S[ and S 2 are robust paths 
for Hi and H2 with parameters rti and ni, respectively. Observe that i < ni + n 2 + 1 must hold. 
Indeed if this would not be the case, then taking any cut with ni + 1 edges in S[ and another cut 
with n 2 + 1 edges in S' 2 yields a cut with ni + n 2 + 2 edges in G, contradicting the fact that S' is 
a robust path with parameter i. We conclude that the algorithm computes S{ correctly in line 23. 
Finally note that the union any two robust paths for the graphs Hi and H 2 with parameters ni and 
n 2 with i < ni + n 2 + 1 yield a feasible solution Si. It follows that the minimum cost such robust 
path is obtained as a minimum cost of a union of two solutions for Hi and H 2 , with robustness 
parameters j and i — j — 1 for some value of j. To allow £3 = Sj or £3 = Sf we let j range from 
— 1 to k and set 5i 1 = = 0. This completes the proof of correctness. 

To prove the bound on the running time, let T(m, k) denote the running time of the algorithm 
on a graph with m edges and robustness parameter k. We assume that the graph is given by a 
hierarchical description, according to its decomposition into single edges. The base case obviously 
takes 0{k) time. Furthermore we assume that the solution (So, • • • , Sk) is stored in a data structure 
for sets, which uses 0(1) time for generating empty sets and for performing union operations. If the 
graph is a series composition then the running time satisfies T(m, k) < T(m', k)+T(m—m' , k)+0(k) 
for some m' < m. If the graph is a parallel composition, then T(m, k) satisfied the same inequality. 
We assume that the data structure, which stores the sets Si also contains the cost of the edges in 
the set. This value can be easily updates in time 0(1) when the assignment into Si is performed. 
It follows that T(m, k) = 0(mk) = 0(nk) as required. □ 

5 Integrality Gap and Approximation Algorithms 

In this section we study the natural fractional relaxation of FTP. We prove a tight bound on the 
integrality gap of this relaxation. This results also suggests a simple approximation algorithm for 
FTP with ratio k+1. We later combine this algorithm with Algorithm Q] to obtain a ^-approximation 
algorithm. 

5.1 Fractional FTP and Integrality Gap 

Let us start by introducing the fractional relaxation of FTP, which we denote by FRAC-FTP. The 
input to FRAC-FTP is identical to the input to FTP. The goal in FRAC-FTP is to find a capacity 
vector x : E — > [0, 1] of minimum cost w(x) = YleeE w e x e sucn t na t for every F € Q(M, k), the 
maximum s-t flow in G — F, capacitated by x is at least one. Note that by the Max-Flow Min-Cut 
Theorem, the latter condition is equivalent to requiring that the minimum s-t cut in G — F has 
capacity of at least one. We will use this fact in the proof of the main theorem in this section. 

Observe that by requiring x € {0, 1} E we obtain FTP, hence FRAC-FTP is indeed a fractional 
relaxation of FTP. 

In the following theorem by 'integrality gap' we mean the maximum ratio between the optimal 
solution value to an FTP instance, and the optimal value of the corresponding FRAC-FTP instance. 

Theorem 12. The integrality gap of FTP is bounded by k + 1. Furthermore, there exists an infinite 
family of instances of FTP with integrality gap arbitrarily close to k + 1. 

Proof. Consider an instance / = (G, M, k) of FTP. Let x* denote an optimal solution to the 
corresponding FRAC-FTP instance, and let OPT = w(x*) be its cost. Define a vector y G M. E as 
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follows. 

f(k + l)x e ife^M 

De = \ (1) 

I minjl, (k + l)x e } otherwise. 

Clearly, it holds that w(y) < (k + 1)0 PT. We claim that every s-t cut in G with capacities y has 
capacity of at least k + 1. Consider any such cut C C E, represented as the set of edges in the cut. 
Let M' = {e € M : x* > ^tpx} denote the set of faulty edges attaining high fractional values in x*. 
Define C = C D M' . If \C'\ > k + 1 we are clearly done. Otherwise, assume \C'\ < k. In this case 
consider the failure scenario F = C . Since x* is a feasible solution it must hold that 

E ( 2 ) 

eeC\C 

Since for every edge e G C \ C" it holds that y e = {k + l)x* we obtain 

E 2/e>^ + l, (3) 
eeC\C 

as desired. From our observations it follows that the maximum flow in G with capacities y is at 
least fc + 1. Finally, consider the minimum cost (k + l)-flow z* in G with capacities defined by 

[ k + l if 6 ^ M (4) 
1 1 otherwise. 

From integrality of c and the minimum-cost flow problem we can assume that z* is integral. Note 
that y e < c e for every e € 22, hence any feasible (/c+l)-flow with capacities y is also a feasible (&+1)- 
flow with capacities c. From the previous observation it holds that w(z*) < w(y) < {k + l)OPT. 
From Lemma [3] we know that z* is a feasible solution to the FTP instance. This concludes the 
proof of the upper bound of k + 1 for the integrality gap. 

To prove the same lower bound we provide an infinite family of instances, containing instances 
with integrality gap arbitrarily close to k + 1. Consider a graph with D ^> k parallel edges with 
unit cost connecting s and t, and let M = E. The optimal solution to FTP on this instance chooses 
any subset of k + 1 edges. At the same time, the optimal solution to FRAC-FTP assigns a capacity 
of -jj^j: to every edge. This solution is feasible, since in every failure scenario, the number of edges 
that survive is at least D — k, hence the maximum s-t flow is at least one. The cost of this solution 
is Taking D to infinity yields instances with integrality gap arbitrarily close to k + 1. □ 

The proof of Theorem [12] implies a simple (k + l)-approximation algorithm for FTP. This 
algorithm simply solves the integer minimum-cost flow problem, defined in proof of the theorem, 
and returns the set of edges corresponding to the optimal integral flow z* as the solution. This 
result is summarized in the following corollary. 

Corollary 13. There is a polynomial (k + 1)- approximation algorithm for FTP. 



5.2 A A;- Approximation Algorithm 

In this section we improve the approximation algorithm from the previous section. The new al- 
gorithm can be seen as a generalization of Algorithm [1] to arbitrary FTP instances. The main 
observation is the following. The reason why the approximation algorithm implied by Theorem [T2l 
gives an approximation ratio of k + 1 is that the capacity of edges in E \ M is set to k + 1, hence, 
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if the flow z* uses such edges to their full capacity, the cost incurred is k + 1 times the cost of 
these edges. This implies that the best possible lower bound on the cost w(z*) is {k + 1)OPTfraCi 
where OPTfrac denotes the optimal solution value of the corresponding FRAC-FTP instance. To 
improve the algorithm we observe that the edges in z* , which carry a flow of k + 1 are cut-edges in 
the obtained solution. 

To conveniently analyze our new algorithm let us consider a certain canonical flow defined by 
minimal feasible solutions. 

Definition 14. Consider an inclusion-wise minimal feasible solution S C E of an instance I = 
(G, s, t, M) of FTP. A flow f s induced by S is any integral s-t (k + l)-flow in G respecting the 
capacity vector 

' 1 if e € S n M 
k + l if e £ S \ M 
if eeE\S. 



To this end consider an optimal solution X* C E to the FTP instance and consider any 
corresponding induced flow /*. Define 

Xpar = {e £ X* : f*(e) < k} and X BRIDGE = {e £ X* : f*(e) = k + l}. 

As we argued before, every edge in Xbridge must be a bridge in H = (V, X*) disconnecting s and 
t. Let u e denote the tail vertex of an edge e £ E. Since every edge e £ Xbridge constitutes an 
s-t cut in H, it follows that the vertices in U = {e u : e £ Xrridge} U {s,t} can be unambigously 
ordered according to the order in which they appear on any s-t path in H, traversed from s to t. 
Let s = «i, • • • , u q = t be this order. Except for s and t, every vertex in U constitutes a cut- vertex 
in H. Divide H into q — 1 subgraphs H , • • • , H q ~ l by letting H l = (V, Y{) contain the union of all 
Ui-Ui + \ paths in H. We observe the following property. 

Proposition 15. For every i £ [q — 1] the set Yi C E is an optimal solution to the FTP instance 
U = (G,Ui,u i+1 ,M). 

Consider some i £ [q — 1] and let /* denote the flow /*, restricted to edges in H l . Note that /* 
can be viewed as a -Uj-iij+i (k + l)-flow. Exactly one of the following cases can occur. Either H l 
contains a single edge e £ E \ M, or 

max/*(e) < k. 

In the former case, the edge e is the shortest Ui-Ui+i path in (V, E \ M). In the latter case we can 
use a slightly updated variant of the algorithm in Corollary [13] to obtain a A;- approximation of the 
optimal FTP solution on instance 1{. Concretely, the algorithm defines the capacity vector 



/ 



k iie^M 
1 otherwise, 



and finds an integral minimum-cost Uj-ttj-fi (k + l)-flow Y* in G, and returns the support Y C E 
of the flow as the solution. The existence of the flow /* guarantees that w(y*) < w(f*), while the 
fact that the maximum capacity in the flow problem is bounded by k gives w(Y) < kw(y*). It 
follows that this algorithm approximates the optimal solution to the FTP instance Ii to within a 
factor k. 

To describe the final algorithm it remains use the blueprint of Algorithm [TJ The only difference 
is in the computation performed in step 5. Instead of finding two edge-disjoint u-v paths, the new 
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algorithm solves the aforementioned flow problem. An appropriate adaptation is also made to step 
12. We summarize the main result of this section in the following theorem. The proof is omitted, 
as it is identical to that of Theorem [61 with the exception of the preceding discussion. 

Theorem 16. There is a polynomial k- approximation algorithm for FTP. 

6 Conclusions and Future Work 

This paper presents FTP, a simple natural extension of the minimum k edge-disjoint path problem, 
interpreted as a fault-tolerant connection problem. In our model not all subsets of k edges can 
be removed from the graph after a solution is chosen, but rather any subset of faulty edges, a 
subset provided as part of the input. Such an adaptation is natural from the point of view of many 
application domains. 

We observed a dramatic increase in the computational complexity of FTP with respect EDP. 
At the same time we identified several classes of instances admitting a polynomial exact algorithm. 
These classes include the case k = 1, directed acyclic graphs and fixed k and series-parallel graphs. 
Next, we defined a fractional counterpart of FTP and proved a tight bound on the corresponding 
integrality gap. This result lead to a ^-approximation algorithm for FTP. We conclude by stating 
a number of promising directions for future research. 

Perhaps the main problem remains to understanding the approximability of FTP. In particular, 
it is interesting to see if the approximation guarantee for FTP can be decreased to the approximation 
guarantees of the best known algorithms for the Steiner Tree problem. It is also interesting to relate 
FTP to more general problems such a Minimum-Cost Fixed-Charge Network Flow and special cases 
thereof. The complexity of fc-FTP in still unknown. It is interesting to see if the methods employed 
in the current paper for 1-FTP and fc-FTP on directed acyclic graphs can be extended to fc-FTP 
on general graphs. It is also natural to study other classical combinatorial optimization problems 
in the model proposed in this paper. In particular, the counterpart of the Minimum Spanning Tree 
problem seems to be a promising next step. We remark that all the results about integrality gaps 
and approximability of FTP can be carried over (with some modifications) to this problem. For 
example, the proof of Theorem [12] can be adapted to the latter problem to give an integrality gap 
of k + 3, in combination with the integrality gap result for the Minimum-Cost fc-Edge Connected 
Spanning Subgraph problem of Gabow et al. |llj . 
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